o 
O 



i 

c 

o 
o 



> 

m 
in 

d 
o 



X 



Power optimization for domain wall motion in ferromagnetic nanowires 

0. A. Tretiakov, 1 Y. Liu, 1 and Ar. Abanov 1 

Department of Physics & Astronomy, Texas A&M University, College Station, Texas 77843-4242, 
USA 

(Dated: October 7, 2010) 

The current mediated domain-wall dynamics in a thin ferromagnetic wire is investigated. We derive the 
effective equations of motion of the domain wall. They are used to study the possibility to optimize the power 
supplied by electric current for the motion of domain walls in a nanowire. We show that a certain resonant 
time-dependent current moving a domain wall can significantly reduce the Joule heating in the wire, and 
thus it can lead to a novel proposal for the most energy efficient memory devices. We discuss how Gilbert 
damping, non-adiabatic spin transfer torque, and the presence of Dzyaloshinskii-Moriya interaction can effect 
this power optimization. 



Introduction. Due to its direct relevance to future 
memory and logic devices, the dynamics of domain walls 
(DW) in magnetic nanowires has become recently a very 
popular topic There are mainly two goals which scien- 
tists try to achieve in this field. One goal is to move the 
domain walls with higher velocity in order to make faster 
memory or computer logic. The other one is inspired by 
the modern trend of energy conservation and concerns a 
power optimization of the domain-wall devices. 

Generally, the domain walls can be manipulated 
whether by a magnetic field^ or electric current^ Al- 
though the latter method is preferred for industrial ap- 
plications due to the difficulty with the application of 
magnetic fields locally to small wires. For this reason, we 
consider in this paper the current induced domain-wall 
dynamics. We make a proposal on how to optimize the 
power for the DW motion by means of reducing the losses 
on Joule heating in ferromagnetic nanowires. 6 Moreover, 
because the averaged over time (often called drift) veloc- 
ity of a DW generally increases with applied current, we 
also address the first goal. Namely, our proposal allows 
to move the DWs with higher current densities without 
burning the wire by the excessive heat and thus archive 
higher drift velocities of DWs. The central idea of this 
proposal is to employ resonant time-dependent current 
to move DWs, where the period of the current pulses is 
related to the periodic motion of DW internal degrees of 
freedom. 

The schematic view of a domain wall in a narrow fer- 
romagnetic wire is shown in Fig. [TJ These DWs are char- 




acterized by their width A which is mainly determined 
by exchange interaction and anisotropy along the wire A. 
Another important quantity is the transverse anisotropy 
across the wire K, which governs the pinning of the trans- 
verse component of the DW magnetization. When no 
current is applied to the wire it leads to two degenerate 
positions of the transverse magnetization component of 
the wall: as shown in Fig. [T] and anti-parallel to it. 

To describe the dynamics of DW in a thin wire we 
derived the effective equations of motion from general- 
ized Landau-Lifshitz-Gilbert^ (LLG) equation with the 
current J, 



S x H, 



J™ 

dz 



oz 



aS x S, 



(1) 



where S is magnetization unit vector, H e ff = 6H/6S is 
the effective magnetic field given by the Hamiltonian H 
of the system, B is non-adiabatic spin torque constant, 
and a is Gilbert damping constant. The derivation of 
the effective equations of motion is based on the fact 
that in thin ferromagnetic wires the static DWs are rigid 
topologically constrained spin-textures. Therefore, for 
not too strong drive, their dynamics can be described 
in terms of only a few collective coordinates associated 
with the DW degrees of freedom^ In very thin wires, 
there are two collective coordinates corresponding to two 
softest modes of the DW motion: the DW position along 
the wire z$ and the magnetization angle <fi m the DW 
around the wire axis. All other degrees of freedom are 
gapped by strong anisotropic energy along the wire. 

By applying the orthogonality condition to LLG, one 
can obtain the equations of motion for the two DW soft- 
est modes, zo(t) and <j>(t)r& 



Zo = A J + B[J - j c sin(2(/>)], 
= C[J-j c sin(2^)], 



(2) 
(3) 



FIG. 1. (color online) A schematic view of a current-driven 
domain wall in a ferromagnetic wire. The DW width is A. 



where J(t) is a time-dependent current. The co- 
efficients A, B, C, and critical current j c can be 

evaluated for a particular model in terms of a, B 

and other microscopic parameters. Following Ref. IToL 
for the model with Dzyaloshinskii-Moriya interaction 
(DMI) one can find A = B/a, B = (a - B)(l + 
oTA)/[a(l + a 2 )], C = (a - B)A/[(1 + a 8 )A§], and 
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FIG. 2. (color online) DW motion characteristics for dc cur- 
rents, (a) Drift velocity Vd of DW as a function of current J 
for B > and B < 0, see Eq. ©. The slope at J < j c is given 
by A, whereas at J > j c it is A + B. (b) Power of Ohmic 
losses Pdc(Vd/V c ) = J 2 /jc as a function of drift velocity Vd- 
For B < the power has a discontinuity at Vd/V c = 1. 



j c = (aKA/ \a - /3\)[nTA/ sinh(7iTA)], where J cx is ex- 
change constant, D is DMI constant, and T = D/J ex . 
Also, A = Ao/Vl - T 2 A^ where A is the DW width in 
the absence of DMI. 

Alternatively, Eqs. © and can be obtained in a 
more general framework by means of symmetry argu- 
ments. We note that because of the translational invari- 
ance Zq and (j) cannot depend on zq. Furthermore, to the 
first order in small transverse anisotropy K, (j) and zq are 
proportional to the first harmonic sin(20). Then the ex- 
pansion in small current J up to a linear in J order gives 
Eqs. ([2]) and (|3j). In this case the coefficients A, B, C, 
and j c have to be determined directly from experimental 
measurements j 1 1 ' 12 

For the dc current applied to the wire the DW dy- 
namics governed by Eqs. ([2]) and ((3|) can be obtained 
explicitly^ For J < j c and A ^ the DW only 
moves along the wire and is tilted on angle <fio from 
the transverse-anisotropy easy axis given by condition 
sin(2</> ) = J/jc- The drift velocity is Vd = (z (J)) = AJ 7 
see Eq. ([2]). Therefore, the linear slope of Vd(J) below j c 
gives constant A, see Fig. [2] (a). The value of j c is deter- 
mined as the endpoint of this linear regime. At J = j c 
the magnetization angle becomes perpendicular to the 
easy axis, 0o = ir/2. For J > j c the DW both moves and 
rotates, and Eqs. © and © give V d = AJ+B-s/J 2 - j%, 
so that the slope of Vd(J) at large J gives A + B. 

Power optimization. The largest losses in the nanowire 
with a DW are the Ohmic losses of the current. In gen- 
eral, the influence of the DW on the resistance is negli- 
gible and therefore we can assume that the resistance of 
the wire is constant with time. Then the time-averaged 
power of Ohmic losses is proportional to (J 2 (t)). Since 
the resistance is almost constant, in this paper we will 
calculate P — (J 2 (t)} and loosely call it the power of 
Ohmic losses. Our goal is to minimize the Ohmic losses 
while keeping the DW moving with a given constant drift 
velocity. 

For the following it will be convenient to introduce the 
dimensionless variables for time, drift velocity, current, 
power, and the ratio of slopes of Vd( J) at large and small 



currents, 

V d . J V A + B 

r = C 3c t, v d = -, j = -, p= ]l , a=—^-. 

(4) 

Although we note that in the special case of a = fi, 
it can be shown that C = B = and one cannot use 
dimensionless variables (|4|). However, in this case the 
DW dynamics is trivial)^ the DW does not rotate 4> = 
0, 7r and moves with the velocity zq = J. 

First, we consider the case of dc current and the power 
as a function of drift velocity. For Vd < 1 we find pdc = 
v\. For currents above j c the power pdc(fd) = j 2 is given 
in terms of drift velocity Vd — j+{B/A)y/ j 2 — 1 as shown 
in Fig.[2](b). The power is quadratic in Vd, and for B < 
it has a discontinuity at v d = 1. 

In general, the DW motion has some period T and 
current j(r) must be a periodic function with the same 
T to minimize the Ohmic losses. Measuring the angle 
from the hard axis instead of easy axis and scaling it 
by 2, i.c, 2<fi = 8 — tt/2, we can write the dimensionless 
current drift velocity as£ 

j(r) =6/2 -cosO, v d = ^{9)-(cos6), (5) 

where 9 = dQ/dr. 

To minimize the power of Ohmic losses we need to find 
the minimum of (j 2 (t)) at fixed Vd, 

p = Ue/2 - cos6>) 2 - 2p{ae/2 -cose- v d )^ , (6) 

where we use a Lagrange multiplier 2p to account for the 
constraint given by Vd from Eq. ([5]). Power ^ can be 
considered as an effective action for a particle in a peri- 
odic potential U , and its minimization gives the equation 
of motion 6/2 = —dU/dd which in turn can be reduced 
to 

6 = ±2y/d-U(6,p), U{0,p) = -cos 2 6 -2pcos6. 

(7) 

where d is an arbitrary constant. Since changing p — > —p 
in U of Eq. ([7]) is equivalent to changing 6 — > it + 6, below 
we can consider only positive p. 

Eq. ([7]) shows that there are two different regimes: 1) 
the bounded regime where d < max[U(#, p)] in which 
case 6 is bounded, and the particle oscillates in potential 
well U(9), see inset of Fig. [3] (a); and 2) the rotational 
regime where d > max[U(#, p)] with freely rotating mag- 
netization in the DW. 

In the bounded regime the particle moves between the 
two turning points — 6q and 9q given by d = U(±6o,p). 
Since 6 is a bounded function (6) = and v d — — (cos 6) . 
One can show-£ that in this regime the power of Ohmic 
losses is minimal for dc current, i.e., p = v% 

In the rotational regime the term in Eq. ([S]) with (9) 
should be kept because 6 is not bounded. The equation 
of motion is the same as for a nonlinear oscillator. 6 Using 
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FIG. 3. (color online) (a) Minimal power of Ohmic losses 
p — (J 2 )/jc as a function of drift velocity Vd shown by solid 
line for a — 0.5. The dashed line depicts p for dc current. The 
inset shows the potential U (8) in which a "particle" is moving 
in the bounded (pendulum-like) and unbounded (rotational) 
regimes. A sketch of (J 2 )(Vd) shown by solid line in (b) for 
/3 » a (a < 1) and (c) for /3 < a (a > 1). 

the minimization condition dp/dp\ Vd — one finds 

/ y / d-U{9,p)d0 = 2-Kap. (8) 

J — 7T 

This equation defines the relationship between d and p. 

The results for the minimal power of Ohmic losses 
p{vd) are presented in Fig. [3l For a > 1 there is a crit- 
ical velocity w rc < 1, such that at Vd < v rc the power 
of Ohmic losses is p — v\ — pdc- Above v IC one can 
minimize the Ohmic losses by moving DW with resonant 
current pulses. Right above i» rc there is a certain range of 
Vd where p = 2poVd — Pq with po(a) < 1 given by Eq. ([8]) 
with d = p 2 . The critical velocity is found as v rc — Pq(cl). 
For a < 1, see e.g. Fig. [3] (a), we find that v IC = 1, 
whereas at Vd > 1 minimal power p is significantly lower 
than pdc- Immediately above Vd = 1 we find that there 
is a range of Vd where p is linear in Vd- At large Vd the 
minimal power is always smaller than p^ c , the difference 
between them then approaches pdc ~P= (1 — l/ a ) 2 /2- 

We note that even in the limiting cases of the systems 
with weak (/3 <C a) or strong (/3 3> a) non-adiabatic spin 
transfer torque, see Fig. [3] (b) and (c), where the power 
of Ohmic losses is high for dc currents, the optimized ac 
current gives dramatic reduction in heating power thus 
greatly expanding the range of materials which can be 
used for spintronic devices We also note that DMI 
suppresses critical current j c and affects parameter a. 

For Vd < v TC the optimal current coincides with the dc 
current, above i> rc the resonant current j(t) is plotted in 
Fig. [4] for a — 2 and two different velocities Vd- At Vd > 
v IC the current's maximum j max increases from 2 — v IC 
at small enough Vd ^ 1 up to j max « Vd/a at Vd 
1. The current's minimum increases monotonically from 
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FIG. 4. (color online) Resonant time-dependent current J(t) 
with t = Cjct for drift velocities Vd = 0.5 (dashed line) and 
Vd ~ 4.5 (solid line) for a — 2. 



small positive values j m j n = v rc at Vd ~ 1 up to j min — 
jmax - 2 1 1 - a\/a at Vd > 1. At u d < 1 (for a > 1) the 
time between the current picks decreases with increasing 
velocity as T ~ (7ra — 2 arcsinw rc — w rc ), whereas the 
pick's width is given by ~ 1-3/a/ (1 — v rc ). Therefore, at 
small Vd — v rc the picks are widely separated, then as Vd 
increases the time between the picks decreases. At Vd ^> 
1 the optimal current has a large constant component 
and small-amplitude ac modulations on top of it. 

Conclusions. We have studied the current driven DW 
dynamics in thin ferromagnetic wires. The ultimate lower 
bound for the Ohmic losses in the wire has been found 
for any DW drift velocity Vd- We have obtained the ex- 
plicit time-dependence of the current which minimizes 
the Ohmic losses. We believe that the use of these res- 
onant current pulses instead of dc current can help to 
dramatically reduce heating of the wire for any Vd- 
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